A multidimensionally consistent version of Hirota's 
discrete KdV equation 



James Atkinson 

School of Mathematics and Statistics, The University of Sydney, NSW 2006, 
Australia. 

Abstract. A multidimensionally consistent generalisation of Hirota's discrete KdV 
equation is proposed, it is a quad equation defined by a polynomial that is quadratic 
in each variable. Soliton solutions and interpretation of the model as superposition 
principle are given, ft is discussed how an important property of the defining 
polynomial, a factorisation of discriminants, appears also in the few other known 
discrete integrable multi-quadratic models. 

PACS numbers: 02.30.1k 
1. Introduction 

This paper mainly concerns the following polynomial that is quadratic in each of four 
variables: 

Q p , q (v, v, v, v) := [p(vv + vv) - q(vv + vv)] 2 

— 4[(v — v) (v — -0) + p — q] [(v — v) (v — v) + (p — q)vvvv] . 

In the simplest setting it defines an equation Q p , q (v , v, v, v) = that relates variables 
on vertices of each quad on the regular Z 2 lattice, whence v = v(n, m), v = v(n + l,m), 
v = v(n,m + 1) and v = v(n + l,m + 1) are values of the dependent variable as a 
function of independent variables n, m e Z, and p,g6C are constant parameters of the 
equation. In the case p — — q — 2 (1) is reducible, 

Q.2,-2(v, v, v, v) = 16(vvv — vvv — v + v)(vvv — vvv — v + v). (2) 

Each polynomial factor in (2) defines a quad equation recognisable as Hirota's discrete 
KdV [1], and it is natural to view (1) as a generalisation of that equation with the 
feature of being multidimensionally consistent. 

The multidimensional consistency identified in [2, 3] is an elegant algebraic 
integrability feature for fully discrete equations. It is a property inherent in systems that 
emerge as the superposition principle for commuting auto-Backlund transformations, 
and therefore models exhibiting the property can be found throughout soliton and 
geometric transformation theory. Most are in the form of rational systems that define 
single-valued evolution from appropriate initial data. A notable exception is a model 
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obtained originally by Kashaev in [4] which relates variables on vertices of a cube in 
Z 3 by a multi- quadratic polynomial, we call it the discrete CKP equation because of 
its interpretation, due to Schief [5], as superposition principle for the CKP hierarchy. 
Probably the only documented example of a quad equation with similar features is the 
discrete KdV equation obtained by Adler and Veselov [6] , 

[{u - u + p + q) (u - u + p - q) - p 2 + q 2 ] 2 

~ ~ ~ ~ ~ ~ ~ \ ) 

— A(u — u) [(u + u){u — u)p — {u + u)(u — u)q + {u — u)pq\ =0, 

which is the superposition principle for Backlund transformations of the KdV equation. 
This model was also derived by Kassotakis and Nieszporski in [7] through an approach 
which allows explicit verification of its multidimensional consistency. As described in 
[6] it is related to the equation 

{w — w)(w — w) = p — q, (4) 

which is superposition principle for Backlund transformations of the potential KdV 
equation [8, 9], by the system 

u + u + p = (w — w) 2 , u + u + q = (w — w) 2 . (5) 

In the continuous setting u = d x w and (5) is the auto-Backlund transformation for the 
potential KdV equation, however in the fully discrete setting the system (5) is viewed 
as a Backlund transformation connecting the distinct equations (3) and (4). The multi- 
affine model (4) is much better studied than its multi-quadratic counterpart (3), this is 
natural because it is simpler and the models are anyway connected due to the Backlund 
transformation (5). Of course it is possible that the class of integrable multi-quadratic 
models, of which (3) is a member, contains systems that are not Backlund related to 
multi-affine ones. For instance such transformation connecting model (1) to a multi- 
afline equation is not presently known. 

The relationship between the equation defined by (1) and the model (3) is actually 
hierarchical, with (3) appearing by the following degeneration procedure 

lim ^—Q 2 ie P ,2ieq(Vi(l + eu), Vi(l + eu), Vi(l + eu), \fi(l + eu)) = 0, (6) 

where i 2 = — 1, as a limiting sub-case of the equation defined by (1). 

The contents of this paper are as follows. The multidimensionality of the new 
example of multi-quadratic model (1) is discussed first in terms of a consistent system 
of polynomials. An important additional property, namely discriminant factorisation, 
is observed in Section 3, it removes the apparent awkwardness of dealing with 
multivalued evolution of the quadratic models and allows straightforward verification 
of the multidimensional consistency. This turns out to be connected to the concept 
of vertex-bond variables discussed by Hietarinta and Viallet in [10]. Soliton solutions 
of the equation defined by (1), which are strikingly similar to the well-known Hirota 
discrete KdV soliton solutions, are given in Section 4. In Section 5 property (2) is used 
to further clarify the relationship between these models. In the final section we discuss 
the importance of the discriminant factorisation property in the context of other known 
multi-quadratic discrete integrable models. 
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Figure 1. Variables assigned to the vertices of a cube. System (8) determines four 
possible values of v from initial data v, v, v and v. 

2. Multidimensionality 

For consistency in multidimensions the key properties of the polynomial (1) are first the 
symmetry 

Qp,q( V i *>» ^> *>) = Qq,p( V i *>» ( 7 ) 

and second the consistency of the system 

Qp,q(v, V, V, V) = 0, Qp,q(v, 5, V, 5) = 0, 

Qq,r(v, V, V, V) = 0, Qg, r (u, V, 5, v) = 0, (8) 
Qr,p(v, V, V, V) = 0, Q r ,p(v, V, V, 5) = 0, 

for any choice of the parameters p, q and r, which is usually visualised by assigning 
variables to vertices of a cube as in Figure 1, and equations to faces. Due to the 
quadratic defining polynomial the consistency needs to be understood by considering 
the three polynomial equations that remain after the elimination of variables v, v and v 
from (8). Direct calculation shows that the remaining polynomials share one irreducible 
factor which is degree four in the remaining variables v, v, v, v and v. This means that 
for generic initial data {v, v, v, v} the system (8) determines four possible values for v. 
An expression for v in terms of this data will be given in the following section. 

These two basic properties, namely covariance and consistency around a cube, imply 
consistency of the d dimensional system 

QiH&( v > T i v > T j v > T i T j v ) = °> e {l---<0, (9) 
where now v is defined on Z rf , Ti . . . are shift operators and pi . . - Pd are corresponding 
lattice parameters associated with each direction. This is the wider and more natural 
setting for the model defined by polynomial (1). Initial data for (9) can be specified on 
co-ordinate axes, but, unlike for the multi-affine defining polynomial, such data does not 
determine a unique solution. However this multi-valued evolution can brought under 
control due to a further important property of the defining polynomial. 
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3. Discriminant factorisation 

It is straightforward to verify that the discriminants of (1) are reducible, specifically 

A[Q p ,q(v, v, v, v),v] = 16 (p — qf{v — v) 2 {l + pvv + v 2 v 2 )(l + qvv + v 2 v 2 ), (10) 

and due to its symmetry there is a similar formula for discriminants of this polynomial 
with respect to the other variables v, v and v. 

The form of (10) suggests introduction of auxiliary variables U\ and cr 2 satisfying 
equations 

a\ = 1 + pvv + v 2 v 2 , <72 = 1 + qvv + v 2 v 2 , (11) 

which are naturally associated with edges of the lattice in the n and m directions 
respectively, see Figure 2. By solving equation Q Pjq (v,v,v,v) = for v and exploiting 
these variables we obtain an equation of the form 

FpJy,V,V,V,(Ti<Td=Q (12) 

which is by construction of polynomial degree one in v and aia 2 , for instance: 

J- p>q {v, v, v, v, 0\0<i) := v\p(vv — vv) + q(vv — vv)\ — 2{v — £>)(ciC2 — 1 — v 2 vv) (13) 

(the precise form of this equation is chosen because it is simple, but it is not unique). 
Sequentially solving for each of the quad variables in the same way leads to the following 

system 

Jp >g (u,U,U,U,<7i<7 2 ) = 0, 

T Ptq (v,v,v,v,<T 1 a 2 ) = 0, ^ 
J r Piq (v,v,v,v,<j 1 <T 2 ) = 0, 

up to the choice of sign of the discriminant terms appearing in the last argument. Each 
equation in (14) implies Q Piq (v,v : v,v) = modulo the relations (11). The choice of 
sign for the discriminant terms is important for self-consistency, ensuring that any one 
equation in (14) is a consequence of the other three. Within this constraint other choices 
of the signs can also break covariance or consistency on the cube, so we restrict attention 
to this canonical one. 

The main benefit of replacing the quad equation defined by (1) with system (14) is 
the single-valued evolution from initial data. The usual initial value problem for a quad- 
equation involves specifying the dependent variable on vertices along some admissible 
lattice path or collection of paths. The modification required for system (14) is that 
variables on path edges, subject to the constraint (11), also need to be specified. For 
instance on the single quad of Figure 2 the variables {v, &i, <7 2 } are determined rationally 
by (14) from {v, v,v, o~i, a 2 }- This system is therefore a model of vertex-bond type as 
described in [10], but with the additional feature of preserving algebraic relations (11) 
on edges. In this setting the multi-valued evolution of the underlying correspondence is 
reflected in the non-uniqueness of initial edge variables when they are defined in terms 
of vertex variables by (11). 
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Figure 2. Variables assigned to the vertices of a quadrilateral, and auxiliary variables 
to the edges. 



As mentioned, the discriminant signs have been chosen in system (14) so that it 
inherits covariance from (1), this enables it to be imposed on each face of the cube in 
Figure 1 without requiring to specify orientation. Initial data to check consistency on 
the cube for (14) is {v, v,v, v, <7i, a 2 , <7 3 } subject of course to the edge relations: (11) 
and a\ = 1 + rvv + v 2 v 2 . Putting system (14) in place of each equation in (8) leads by 
direct calculation to the following equation: 

v(vv — vv)r(p — q) 2 — 2(v — v)(aicr 2 — 1 — v 2 vv)(r — p)(q — r) 

+v(vv — vv)p(q — r) 2 — 2(v — v)(a 2 a 3 — 1 — v 2 vv)(p — q)(r — p) (15) 

+v(i)v — vv)q(r — p) 2 — 2(t> — v)(a 3 ai — 1 — v 2 vv)(q — r)(p — q) = 0. 

This determines v uniquely in terms of the initial data on the cube. Transformations 
(<7i, a 2 , 03) (±o"i, ±o" 2? ±o"3) yield other polynomial roots of the related multi- valued 
system (8). Note that although this transformation group contains eight elements, 
consistent with our earlier claim it only leads to four different values of v due to the 
symmetry (a u cr 2 , cr 3 ) ->■ {-a u -cr 2 , -cr 3 ) of (15). 

4. Multi-soliton solutions 

The consistency on the cube provides immediately an auto-Backlund transformation for 
the equation defined by (1): 

Qp,k(v, v, v, v) = 0, Qq,k(y, v, v, v) = 0. (16) 

This system determines a new solution v = v(n,m) from old solution v = v(n,m), 
k is the Backlund parameter. The Backlund equations determining v are biquadratic 
correspondences associated with edges of the lattice. 

That the equation defined by (1) admits the constant solution v — c can be seen by 
inspection, and it is natural to ask for the multi-soliton solution obtained by subsequent 
iteration of the Backlund transformation. At the first step of the Backlund chain 
the equations (16) are compatible autonomous symmetric biquadratic correspondences 
whose solution is standard. Superposition can then be used to obtain the next few 
solutions. Comparison with solutions obtained in [1] motivates a slight generalisation 
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of Hirota's original substitution and leads to the following 7V-soliton solution for model 

v = cf FT- (17) 

Here / and its counterparts are standard Hirota-type polynomials 

f : = E ^11^ f + ' : = E YjJlatarfH, 

ic{i...n} iei ic{i...n} iei ^ 

f +:= E Y iH a tf>i> f~ : = E Y iU a iP- 

IC{1...N} iei IC{1...N} iei 

where constants Yj, I C {1 . . . N} and af . . . are defined as 

Y * - 11 A^T^J ' ° l -lic^ + fc?' (19) 

functions p\. . . are defined as 



,P* + k*J \q* + k* 

and parameters p*, g* and k\ . . .k* N are related to the lattice and Backlund parameters 
respectively by equations 

pl = l +P c 2 + c\ ql = l + q c 2 + c\ {klf = 1 + k t c 2 + c\ ie{l...N}. (21) 

Constants pi,o ■ ■ ■ Pn,o an d parameters k\...k* N can be chosen freely in the solution. 
Only the solution is given here, we omit recalling standard methods required to verify 
it (cf. [11, 12, 13]). By inspecting (18), (19), (20) and (21) above it can be seen that 
the specialisation {p, q} = {+2, —2} returns solution (17) to the one obtained originally 
by Hirota: 

v = c££, (22) 
// 

which is consistent with observation (2) made earlier. 



5. Relationship with Hirota's discrete KdV equation 

In this section model (1) will be used to obtain a Backlund transformation for the 
equation due to Hirota [1]. The method relies on the multidimensional consistency of 
the associated single- valued vertex-bond system (11), (14), combined with property (2). 
Note that the transformation obtained by this method should be identified with the 
bilinear Backlund transformation given in [1]. 

The specialisation p = — q = 2 reduces the system (11), (14) to 

<7i = (-l) n+m (l + w), a 2 = (-l) n+m (l - vv), (23) 
-11 
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up to transformations (o"i,o- 2 ) — > (±<Ti, ±0-2), in particular recovering Hirota's equation 
(24). With the same specialisation the Backhand transformation for (11), (14), which 
is inherent from its three-dimensional consistency, leads to the following Backlund 
transformation for (24): 

v(vv — vv)(r — 2) = 2(v — v)(l + vv)(a 3 — 1 — vv), 

v(vv — vv)(r — 2) = 2{y — v)(l + vv )(a 3 + 1 + vv), 

v(vv — vv)(r + 2) = 2{v — v)(l — vv)(a 3 — 1 + vv), 

_ ± . (26) 

v(vv — vv)(r + 2) = 2(v — v)(l — vv)(a 3 + 1 — vv), 

where substitution a 3 — > (— l) n+ "V 3 has also been made in order to obtain the 
transformation in autonomous form. 

For a fixed function v = v(n,m) the systems (25) and (26) define bi-rational 
mappings (^,03) h-> (v,cr 3 ) and (^,03) h-> (w, <t 3 ) which preserve the relation a\ = 
1 + rvv + f 2, u 2 , and which are associated with edges of the lattice in the n and m 
directions respectively. The constraint (24) on v is sufficient for compatibility of the 
mappings, and they determine a new function v = v{n, m) that again satisfies (24). The 
free parameter r is the Backlund parameter. 

This Backlund transformation determines not only v but also the auxiliary variable 
(T3, which then plays a role in the superposition principle. The superposition principle 
itself is obtained from the four-dimensional consistency of the system (11), (14): let 
(v, <t 3 ) and (v, <r 4 ) be obtained by the Backlund transformation (25), (26) from a solution 
v of (24) with the corresponding Backlund parameters being r and s. Then the equation 

F r , a {v, v, v, v, a 3 (j 4 ) = (27) 

determines a new solution v = v(n, m) of (24) related to v by the composition of the two 
transformations. By using edge relations a\ — 1 + rvv + v 2 v 2 and a\ — 1 + svv + v 2 v 2 
the superposition principle (27) can be expressed without the auxiliary variables as 
Qr, s (v, v, v, v) = 0. Thus the full parameter-dependent model (1) is recovered as 
superposition principle for the Backlund transformation (25), (26) of Hirota's equation 
(24). 

6. Discussion: multi-quadratic models 

To conclude we discuss further the discriminant factorisation property and method of 
introducing auxiliary variables which was used in Section 3. 

Weierstrass biquadratic 

The motivating example for the technique is classical theory of the Weierstrass 
biquadratic equation: 

(uu + up + pu - eie 2 - e 2 e 3 - e 3 e ± ) 2 

(28) 

—A(u + u + p - ei - e 2 - e 3 )(uup - t\e 2 e 3 ) = 0. 
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The single-valued system equivalent to the quadratic model (28) is usually written as 
u + u + p — e 1 — e 2 — e 3 , (u — p)U + (u — u)P — {p — u)U, (29) 



1 {U -P^ 2 



4 y u — p 

where the auxiliary variable and parameter satisfy 

U 2 = 4(u - ei )(u - e 2 )(« - e 3 ), P 2 = 4(p - ei)(p - e 2 )(p - e 3 ). (30) 

System (29) defines a single-valued mapping (u, U) i-> (w, C/) which preserves the 
algebraic relation satisfied by the auxiliary variable in (30). In this setting the 
analogue of multidimensional consistency is the well-known commutativity between the 
correspondence (28), or the associated single- valued mapping (29), and its counterpart 
with different choice of parameter p. 

Adler-Veselov model and the Kassotakis-Niesporski Idea-system variables 

The multi-quadratic model (3) was re-discovered in [7] through an associated single- 
valued Idea system involving variables only on lattice edges. This led in particular to 
a single- valued vertex-bond system equivalent to (3), and we wish to briefly show how 
the same system can be recovered by the constructive method based on discriminant 
factorisation. 

The discriminant, with respect to variable u, of polynomial equation (3): 

16(p — q) 2 {u — u) 2 {u + u + p){u + u + q) , (31) 
motivates introduction of the edge variables 

o\ = u + u + p, o\ = u + u + q, (32) 

which are exactly the variables of the associated Idea system. As in Section 3 a single- 
valued vertex-bond system replaces (3) after introduction of the edge variables (32), in 
this case it can be written as follows: 

(p + q)(u — u) — (p — q){u — u) — (u — u){2aiU2 — 2m — u — u), 
(p + q)(u — u) — (p — q){u — u) — (u — u)(2aia2 — 2u — u — u), (33) 
(p + q){u — u) — (p — q){u — u) — (u — u)(2aicr2 — 2u — u — u). 
Note that the Idea system itself, which may be written in the form 

o~i + a 1 = a 2 + cr 2 , 0i<7i —p — o 2 o 2 - q, (34) 

is also a consequence of system (33) and edge relations (32). Direct calculation using 
the vertex-bond system (33) was used in [7] to verify consistency on the cube of (3), we 
give here the same formula in a slightly different form more similar to (15): 

(u + u — u — u)r{p — q) 2 — (u — u)(2aia 2 — 2u — u — u){r — p)(q — r) 

+(u + u — u — u)p(q — r) 2 — (u — u)(2a2cr 3 — 2u — u — u)(p — q)(r — p) (35) 

+(u + u — u — u)q{r — p) 2 — (u — u)(2a 3 ai — 2u — u — u)(q — r){p — q) =0, 

where o~\ — u + u + r. Note that the symmetry (ci, cr 2 , cr 3 ) — > (—o~\, — cr 2 , — o" 3 ) of (35) 
means that initial data {u, u, u, u} determines only four distinct values of u. 
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The discrete CKP equation 

The symmetric discrete Darboux system was given by Schief in [5] as a single-valued 
system involving variables on quad faces of the cubic lattice Z 3 (cf. also [14, 15]). 
Introduction of an associated function r on vertices as a potential for the face variables 
led to the following equation defined by a multi-quadratic polynomial: 

(tt — ff + ff — ff) 2 — 4(tt — ff)(ff — ff) = 0, (36) 

which was shown to be superposition principle for the CKP hierarchy, and therefore it 
should be consistent in multidimensions. This was confirmed in calculations by Tsarev 
and Wolf in [16]. Here we observe that a single- valued system can be constructed from 
(36) due to it having the discriminant factorisation property. The discriminant, with 
respect to f , of the multi-quadratic defining polynomial in (36) is simply 

16(ff — tt) (ff — tt)(tt — tt ), (37) 

and similar expressions hold for discriminants with respect to other variables due to 
(slightly hidden) symmetry of this polynomial. The discriminants motivate introduction 
of new variables, this time on quad faces: 

°"l2 = TT - TT, 0"23 = TT ~ TT i °~13 = ff - TT. (38) 

The equivalent system with single-valued evolution that emerges is: 

-(TT - TT - TT - TT) + (712(723(713 + TTT = 0, 

-(tt -tt -tt - ff) - a 12 a 23 a 13 + rff = 0, 

f , - m) 

-(fT -TT -TT - fr) - (712(7230-13 + TTT = 0, 

T * ~ - ~ ~- 

-(ff - ff - fT - Tf) - Oi 2 (723013 + TTT = 0. 

It contains mixed face and vertex variables unlike the original symmetric discrete 
Darboux formulation of this model which involves variables only on faces. In this 
respect this formulation is not so natural, however it demonstrates again the constructive 
method based on the discriminant factorisation property. And for instance it allows 
direct verification of the multidimensional consistency of (36) as follows. 

First it is easily verified that (36) is covariant, i.e. invariant under interchange of 
lattice directions, and this is inherited by (39). It therefore remains to verify consistency 
of system (39) on (cubic) faces of a four-dimensional hypercube, it requires defining 
additional auxiliary variables 

(7 14 = TT - TT, (T 24 = TT - TT, a 34 = TT - TT, (40) 

and a new shift (grave) in a fourth lattice direction. Direct calculation obtains the 
equation 

t\ I ■ ^ -x, 

— (TT - TT - TT - TT) + (J 12 (T 2 ? y C> 1Z T + O 23CT 34a ' 24 T + (7i 3 (Ji 4 (73 4 T + ^^Cl^ 
+0"120"340'130"24 + 0" 13 Cr 2 4(7l40"23 + ^14^23^12^34 + ffff = 
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which determines r symmetrically from initial data on the hypercube. Note that 
consistent with the observation made in [16], initial data {r, f, f, f , f, f, f , f, f, f, r} 

determines only eight values of f. This can be seen here because although the 
transformation group (a 12 , cr 13 , a u , cr 23 , cr 24 , cr 34 ) ->• (±cr 12 , ±<7i 3 , ±cr 14 , ±cr 23 , ±cr 24 , ±cr 34 ) 
contains 64 elements, the equation (41) has four symmetries of the form 
(o"i2, 0-13, o-i4, o"23, ^24, 034) (-o~i2, -^13, -<7i 4 , cr 23 , ct 24 , ct 34 ) which generate a sub-group 
of eight elements. 

Conclusion 

It has been observed here how a discriminant factorisation property can assuage the 
apparent awkwardness associated with the multi-valued evolution of multi-quadratic 
discrete models. Moreover, its ubiquity amongst the known integrable examples suggests 
to investigate the extent to which this simple feature is connected to their integrability. 
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